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ABSTRACT 

We present a new method aimed to handle long Gamma-Ray Burst (GRBs) as cosmic rulers. 
The recent discovery of a tight correlation between the collimation corrected GRB energy 
and the peak of the 7-ray spectrum has opened the possibility to use GRBs as a new cate- 
gory of standard candles. Unfortunately, because of the lack of low-z GRBs, up to now this 
correlation is obtained from high-z GRBs with the consequence that it depends on the cosmo- 
logical parameters we pretend to constrain. Hopefully this circularity problem will be solved 
when, in a few years, the low-z GRB sample will be increased enough. In the meanwhile we 
present here a new Bayesian method that eases the aforesaid circularity problem, and allows 
to introduce new constrains on the cosmological (JIm, ^^a) diagram as well as to explore the 
universe kinematics up to z w 3. The method we propose offers the further advantage to make 
handy the problem of the (JIm, ^a) loitering line singularity which inevitably appears when 
standard candles with z > 2 are used. The combination of GRB with SN la data makes the 
popular ACDM cosmology more consistent with the Hubble diagram at a 68% confidence 
level. For a flat cosmology we find f^M = 0.28 ± 0.03 for the combined GRB+SN la data 
set. Correspondingly, the transition redshift between cosmic deceleration and acceleration is 
2;t = 0.73 ±0.09, slightly larger than the value found by considering SNe la alone. We briefly 
discuss our results also in terms of non-ACDM dark energy models. 

Key words: cosmological parameters — cosmology :observations — distance scale — gamma 
rays: bursts 



1 INTRODUCTION 

The Hubble diagram of type la Supemovae (SNe la) has revolu- 
tionized cosmology, showing that the current expansion of the uni- 
verse is accelerated (Riess et al. 1998; Riess et al. 2004; Perlmutter 
et al. 1999), with a possible transition from deceleration to accel- 
eration at low redshifts (Riess et al. 2004; Dicus & Repko 2004). 
The most popular explanation for the recent acceleration implies 
the existence of a dominant energy density component with nega- 
tive pressure (dark energy, see for a recent review Sahni 2004). Its 
physical nature is still a mystery. The simplest candidate for dark 
energy is the modern version of the Einstein cosmological con- 
stant A, which, however, raises several theoretical difficulties (see 
Padmanabhan 2003; Peebles & Ratra 2003; Sahni 2004). There is 
now a flurry of activity focused to unveil the nature of dark energy 
or to constrain any alternative explanation to accelerated expan- 
sion. The forthcoming SNAP mission ^ (^max ~ 1-7) will offer 
an unprecedented accuracy on this undertaken. Besides accuracy, 
also a survey reaching larger redshifts is desirable in order to over- 
come parameter degeneracies, and to study the evolution of dark 
energy (Linder & Huterer 2003; Weller & Albrecht 2002; Nesseris 

* E-mail: firmani@merate.mi.astro.it 
^ http://snap.lbl.gov/ 



& Perivolaropoulos 2004). Furthermore, alternative standard can- 
dles to SNe la, free of extinction issues, are also highly desired. 
Long GRBs have been pointed out as possible candidates to be 
these high-redshift candles (Schaefer 2003), and in Ghirlanda et 
al. (2004b, GOLF thereafter) we have shown that they are indeed. 

Assuming isotropy, the energy emitted by GRBs in 7-ray s 
(Eiso) is spread over 4 orders of magnitude, but the presence of an 
achromatic break in their afterglow lightcurve is a strong evidence 
that their emission is coUimated into a cone of semiaperture angle 6 
(e.g. Rhoads 1997; Sari, Piran & Halpem 1999). Correcting for this 
anisotropy, the resulting 7-ray energetics E-, — (1 — cos^)£'iso 
clusters around E-y ~ 10^^ erg, with a small dispersion (0.5 dex. 
Bloom, Frail & Kulkami 2003; Frail et al. 2001), yet not small 
enough for a cosmological use (Bloom, Frail & Kulkami 2003). 

Recently, a very tight correlation has been found (Ghirlanda, 
Ghisellini & Lazzati 2004a, GGL thereafter) between E-y and the 
peak energy -Bpcak of the prompt emission (in a i/F^. representa- 
tion) (§2). This correlation has been used for a reliable estimate of 
-Eiso, and therefore of the luminosity distance, with the aim to con- 
struct a Hubble diagram up to z = 3.2 (GGLF; Dai, Liang & Xu 
2004). 

The existence of the E^-Ep^ak correlation, allowing to know 
the intrinsic emission of GRBs, makes them standard candles sim- 
ilar to SNe la and Cepheid stars. Hopefully, this correlation will be 



2 Firmani et al. 




Figure 1. Rest frame correlations between the spectral peak energy -Epcak and the collimation corrected energy E~i (left panel), the collimation angle 
6 (middle) and the isotropic bolometric energy (right), for the 15 GRBs considered in GGL, assuming nM=0-3 and f!A=0-7- Different symbols (colors) 
correspond to different redshift ranges: circles (blue), triangles (green) and squares (red) are for z < 0.7, 0.7 < z < 1.55 and z > 1.55, respectively. The 
solid line in the left panel represents the best fit i?pcak^^7 correlation (x^ = 1.27) and the shaded region its uncertainty computed in the coordinate frame 
where the variances on its parameters are uncoiTelated. We also show (open circle) the recent GRB 041006 (data in Galassi et al. 2004, D'Avanzo et al. 2004 
and the HETE2 web page i http://space.mit.edu/HETE/Bursts/GR B041006l . 



calibrated from low-2: GRBs, and corroborated by a robust theo- 
retical interpretation. Currently low-z GRBs are scarce: we must 
wait for a few years before having a sufficient number of GRBs at 
z <0.2. For the time being, the correlation is obtained from high- 
2 GRBs with the consequence that it depends on the cosmologi- 
cal parameters which we pretend to find. We present here a new 
method with respect to GGLF that eases this difficulty making pos- 
sible the full use of the available information. We will assume that 
the ACDM cosmology has h = 0.71, = 0.27, and = 0.73. 



2 THE E^-EpEAK CORRELATION 

In Fig. 1 (left panel) we show the correlation between the rest frame 
^Jpeak and Ej for the 15 bursts considered by GGL, using the data 
listed in their Tab. 1 and Tab. 2, and assuming Q,m ~ 0.3 and 
SIa = 0.7. With this choice of cosmology, the reduced ~ 1-27 
for the fit with a powerlaw. The mid panel shows iJpcak as a func- 
tion of the jet aperture angle 6, while the right panel shows Bpeak 
as a function of the isotropic emitted energy Eiso (the "Amati" re- 
lation, Amati et al. 2002) for the same 15 bursts. Also shown is 
the recent GRB 041006: although fully consistent with the GGL 
correlation, it is not included in the sample of 15 GRBs used for 
cosmology due to the still unpublished uncertainties on the relevant 
parameters. The shaded area in the left panel shows the region of 
uncertainties of the best fit correlation: i5peak/(263 ± 15 keV) = 
[E.,/{4:.2 X 10^''erg)]° ™±° °*^ The errors on its slope and nor- 
malization are calculated in the "baricenter" of Speak and E^, 
where the slope and normalization errors are uncorrected (Press 
et al. 1999). For similar i?pcak, there is a distribution of jet aperture 
angles, which is why E.^, for a given -Epcak, is much less spread 
than the corresponding Eiso values. 

In the present sample of bursts with known redshift 
GRB 980425 (Galama et al. 1998) and GRB 031203 (Gotz 
et al. 2003) are outliers for both the Speak — Eiso and the 
E^ — -Epcak correlation. These events are indeed peculiar for their 



small prompt and afterglow emission, as recently discussed by 
Soderberg et al. (2004a). It has been also suggested that these two 
events are sub-energetic, with respect to classical GRBs, because 
they are observed off-axis (Waxman 2004; Yamazaki, Yonetoku & 
Nakamura 2003). Also GRB 020903 (Soderberg et al. 2004b) may 
be an outlier for the E.y — Epcak correlation, although, as discussed 
in GGL, its afterglow data are not conclusive. 

3 THE METHOD 

A Friedmann-Robertson- Walker cosmological model is defined by 
the Hubble parameter H{z) given by: 

H\z) = [nuil + zf + nfe(l + zf + QEfiz)] (1) 

where = 1 — — ^e, J^m and JIe are the present-day den- 
sity parameters of matter and dark energy components, and f{z) 
is a function that depends on the equation-of-state index of the 
dark energy component, w{z) — pe{z) / p-E{z)c'^ . Later on we 
will take into account the case of a cosmological constant (A- 
models) with JIe = ^a, w = —1 and f{z) = 1, as well 
as a linear parametrization of the dark energy equation of state 
(Wi-models) adopting w{z) — wo + w'z, a flat cosmology, and 
f{z) = (l + 2)^'^+'"o-"'')exp(3«;'z). Once defined //(z), the lu- 
minosity distance as a function of z (Hubble diagram) is calculated 
straightforwardly. The transition redshift zt from deceleration to 
acceleration for A-models is given by 2flA = f2M(l + zt)^, and 
for Wi-models by Qm = (1 — nM)( — 1 — 3wo — 3w'2t)(1 + 
2^)3(»o-»')exp(3i(;'2T). In this letter we shall fit observations 
with only two free parameters. These ones are (i) (JlMif^A) for 
A-models, representative of the cosmological constant models, (ii) 
(JIm , Wo) with a flat cosmology and w' — for Wa-models, which 
identify quiessence dark energy models, and finally (iii) (wo,w') 
with (r^M = 0.27, ^A = 0.73) for Wi-models, to identify ki- 
nessence models. 
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The use of GRBs as cosmic rulers brings up two major diffi- 
culties: 

1) Due to the current poor information on low-jz GRBs, the 

— -Epcak correlation (and its scatter) depends on the assumed 
cosmology. This difficulty hides a circularity problem and has 
to be handled with care otherwise it leads to a tautology: if a 
cosmology is used to fix the — iSpcak correlation then the 

— -Epcak correlation cannot be used to find a cosmology. This 
problem was not taken into account by Dai, Liang & Xu (2004). 

2) Because of the previous point and the complex behavior of 
the luminosity distance dL as a function of z, JIm, J^e, some math- 
ematical undesirable attractors appear. This effect can be seen in 
Fig. 1 of GGLF where the GRB confidence level curves are tan- 
gent near S1m=0.07 and r2A=1.15. The effect is basically due to the 
rather singular behavior of near the loitering line that discrimi- 
nates between Big-Bang and no Big-Bang universes (close to the 
region probed by GRBs, but not by SNe). Here cLl falls to zero for 
2 > 1 . The sharpness of this fall increases with z. The obvious result 
is that for high z the constant dh lines tend to be tangent to the loi- 
tering line converging towards S1m~0 and S1a~1. Due to the high 
GRB redshifts, this is just the region where GRBs lie (this problem 
is practically irrelevant for SNe la). If the statistical approach is not 
sufficiently powerful (as the method used in GGLF), the singu- 
larity of is able to attract the most probable solution to some 
improper value (for the method used in GGLF this solution is, 
in fact, f^M^O.l and JIa^LI). 

It is necessary to introduce a new powerful statistical approach 
in order to avoid mathematical artificial tricks in the GRB case. The 
Bayesan approach we propose below solves the difficulties stated 
before and it might come in useful in other similar problems. 

In the case of A-models, for any given cosmology Sl = 
{CIm, ^a) there is a "best fitted" E^-Epcak correlation. For each 
GRB this correlation allows to calculate an observed dLi{0,) = 
[_Biso(l + z)/{4nF)]^/'^, where F is the fluence and Eiso is de- 
rived following Tab. 1 , Tab. 2 and Eq. 1 of GGL. Now, on the Hub- 
ble diagram, for each cosmology we can define X^i^,^) = 
'^i-iid-Lii^) — dhizi, Q))^ /af, where dL{zi, fi) is the theoret- 
ical luminosity distance, cTi the standard deviation on dLi(n), and 
the sum is extended on the A'^ GRBs of the sample. Making use 
of the incomplete gamma function (Press et al. 1999), we trans- 
form -x^{Q,,Cl) into the conditional probability P(r2|n) that pro- 
vides the probability for each Q, given a possible A-defined corre- 
lation. In order to understand the essential of the method, we can 
imagine that at first all Sis (i.e., the corresponding correlations) are 
assumed to be equally probable. Thanks to the previous conditional 
probability, it follows that we do have some knowledge of which is 
the probability of Q,. This is derived from the same P(f2|n) adding 
the contributions of all possible H-defined correlations. Taking this 
probability we assign now a weight to f2, then we repeat the cycle. 
We arrive to another set of probabilities for Q.. We continue to it- 
erate, until convergence. This intuitive approach find a simple and 
elegant formalization through a Bayesian approach by the formula 

p{n) = J p{n\Ci)p'{Ci)dCi (2) 

where the integration is extended on the accessible Q, plane. The 
prior probability P' defines the weights for the £7 -Epcak correla- 
tion, the likelihood P{^l\Cl) is the conditional probability defined 
above, and P is the posterior probability that measures the good- 
ness of fit of the possible cosmologies. We convert the previous 
formula into an equation setting P' — P. A Montecarlo method 
provides an economical approach to find the numerical solution of 
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Figure 2. Constraints in the S^M-f^A plane derived from our GRB sample 
(15 objects, red contours), from the "Gold" SN la sample (156 objects; blue 
solid lines, derived assuming a fixed value of Ho = 65 km s~^ Mpc^^, 
making these contours slightly different from Fig. 8 of (Riess et al. 2004), 
and from the subset of SNe la at 2 > 0.9 (14 objects, blue long dashed 
lines). The three colored ellipses are the confidence regions (orange: 68%; 
light green: 90%; yellow: 99%) for the combined fit of type la SN+GRB 
samples. Dashed lines correspond to the changing sign of the cosmic ac- 
celeration [i.e. q{z) = 0] at different redshifts, as labelled. Crosses are 
the centers of the corresponding contours (red: GRBs; blue: SNe la, white: 
GRB+SN la). The black dot marks the ACDM cosmology. The dotted line 
corresponds to the statefinder r = 1, in this case it coincides with the flat- 
ness condition. 

this equation. For M^o-models the previous procedure is applied to 
the parameter pair {Q,m,'wo), and for M^i -models the same is ap- 
plied to (mo, w'). 

This approach is quite different from mapping the param- 
eter for all points in the fl plane by minimizing the scatter of the 
data points around a f2-dependent correlation (GGLF, see below). 



4 RESULTS 

For fitting the cosmological parameters we have selected the 15 
GRBs reported on GGL (see §2). The GRB sample extends from 
z — 0.16 to 2 = 3.2. The SN la sample, that we combine with 
the GRB sample, includes the most recent discovered high redshift 
SNe la (up to 1.7) and is composed by 156 objects (the "Gold" 
sample, Riess et al. 2004). 

For the A-model, we constrain Qm and SIa. Fig. 2 shows the 
confidence levels (CL) for (i) GRBs alone (adopting our Bayesian 
method), (ii) SNe la alone (for the entire "Gold" sample and for 
only the subset of 14 SNe la at 2 > 0.9, using in both cases the 
standard approach), and (iii) the combination of all SNe la and 
GRBs. The dotted line shows the statefinder parameter r — 1 
(Sahni et al. 2003; see also the jerk parameter in Visser 2004, and 
Riess et al. 2004) that in A-model coincides with the flatness con- 
dition. 
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Figure 3. Constraints on and Hm for a flat cosmology with darlc energy 
wliose equation of state is constant (wq = — 1 corresponds to a cosmolog- 
ical constant). To construct tlie GRB contours we liave followed the same 
method as in Fig. 2. Colored regions: combined (GRB+SN la) constraints 
(color code as in Fig. 2). Dashed lines correspond to the changing sign of 
the cosmic acceleration [i.e. q{z) = 0] at different redshifts, as labelled. 
Crosses are the centers of the SNe la contours (blue: SNe la alone, white: 
GRB+SN la). The black dot marks the ACDM cosmology. The dotted line 
corresponds to the statefinder r = 1. 



The difference between Fig. 2 and Fig. 1 of GGLF is remark- 
able. In Fig. 2 the most probable solution is shown by the red cross 
which corresponds to Q.m=Q-^1 and r2A=0.91, while in Fig. 1 of 
GGLF the best fit model corresponds to f2M— 0.1 and JIa — l-l- 

For GRBs the CL is a little larger than for high-z SNe la, 
but this CL is expected to be reduced by new experiments. Note 
that the orientations of the CL contours are related to the average 
redshift of the sample; the increase of this redshift from all SNe la 
to high-2: SNe la alone, and from high-z: SNe la to GRBs produces 
a counterclockwise rotation of the CL contour pattern. 

When we combine SNe la with GRBs the best-fit model 
moves from Q.m = 0.48, Q,h = 1-00 for the SN la sample to 
r^M = 0.37, r^A = 0.87 for the combined SN+GRB sample, 
favouring a closed universe. However, due to the high correlation 
of the CL contours in Fig. 2, the best-fit value does not mean 
much. Instead, it is worth to remark that the flat geometry solution 
JIm +nA = 1 is compatible now with the 68% CL. Assuming a flat 
geometry prior we obtain Q,m = 0.28 ± 0.03, in excellent agree- 
ment with independent galaxy clustering measurements (Hawkins 
et al. 2003; Schuecker et al. 2003) and the ACDM cosmology. In 
Fig. 2 we also plot the curves of constant zt. For our best-fit model, 
zt = 0.67, while for the flat geometry case our constraints give 
ZT = 0.73 ± 0.09. 

For Wo-models Fig. 3 shows the constraints on wo vs Q,m (i-e. 
a constant equation of state). Best values for the SN-l-GRB sam- 
ple are f^M = 0.44 and wq = —1.68 with zt = 0.40, while 
the ACDM is compatible with the 68% CL (it was not for SN 
la alone). Inside the uncertainty, wq < —1 {r > 1) quiessence 
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Figure 4. Constraints on the wq and w' parameters entering the dark energy 
equation of state (see §3; idq = — 1 and w' = correspond to the cos- 
mological constant Qa)- We assume a flat geometry and Vlyi=0.21. Blue 
contours: constraints from type la SNe; red contours: constraints from our 
GRBs. To construct the GRB contours we have followed the same method 
as in Fig. 2. Colored regions: combined (GRB+SN la) constraints (color 
code as in Fig. 2). Dashed lines correspond to the changing sign of the cos- 
mic acceleration [i.e. q(z) = 0] at different redshifts, as labelled. Crosses 
are the centers of the SNe la contours (blue: SNe la alone, white: GRB+SN 
la). The black dot marks the ACDM cosmology. Dotted line corresponds to 
the statefinder r = 1. 

dark energy models appear favoured. However, while SN la best- 
fit model has uiq < —1, GRBs alone prefer too > —1, even if 
close to —1. A more complete sample of both objects is expected 
to support ACDM. For VFi -models Fig. 4 shows the constraints on 
the dark energy parameters wq vs w' . For the SN+GBR sample the 
best-fit values are too = —1.19 and to' = 0.98 with zt = 0.55. 
The ACDM again is now fully within the 68% CL of the com- 
bined SN+GRB fit. Similar to the previous case even now, inside 
the uncertainty, r > 1 kinessence models (i.e. Chaplygin gas, type 
1 braneworld models, see Alam et al. 2003) seem favoured. In both 
Wi, SN la alone marginally reject ACDM cosmology. This ten- 
dency is balanced by GRBs making hopeful a ACDM cosmology 
for the combined sample. We have here a clear example how GRBs 
complement the cosmological information derived by SN la. 



5 DISCUSSION 

Several kinematical and dynamical parametrizations of cosmology, 
aimed to optimize the fittings to the available observations, have 
been proposed. However current Hubble diagrams of SNe la need 
to be extended to higher redshifts to constrain such parameters with 
better accuracy. 

GRBs are the natural objects able to extend the measure of the 
universe to very high redshifts. This is now possible thanks to the 
Speak vs. E-y empirical correlation, that allows to estimate Eiso of 
each GRB. At present this correlation is known on the ground of 15 
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Figure 5. Residuals of the distance moduli of SNe la and of our GRBs 
with respect to the ACDM case f!M=0.27, nA=0.73. Also shown are 
the differences of various other cosmological models with respect to the 
Om=0.27, f2A=0.73 case: red solid line: Qm=0-2, f2A=0-8; orange solid 
line: f2M=0.37, nA=0.85; blue dot^ashed line: r2M=0.45, Qa=0.95; 
black dashed line: f2M=l, f^A=0; black dotted hne: Om=0.01, S1a=0.99. 
The arrow marks the redshift of GRB 000131 (2=4.5) which is the most dis- 
tant GRB of known z. Since for this GRB there is no measured jet-break 
time, it is not included in our sample. 



GRBs distributed on redshifts up to 3.2; therefore its knowledge de- 
pends on the assumed cosmology. We presented here a new method 
based on a Bayesian approach that takes into account our lack of 
knowledge on the correlation, and optimizes the information to de- 
rive the best-fit model and the CL contours. 

Despite the small number of useful GRB events and their 
somewhat large observational errors, a clear trend emerges when 
combining GRB and SNe la data. The inclusion of GRBs makes the 
ACDM cosmology compatible with the 68% CL. A similar trend, 
but less pronounced, is observed for the SN la sample alone when 
including and not the 16 HST high-2: SNe (see also Choudhury & 
Padmanabhan 2004; Alam et al. 2004). 

Fig. 5 shows the Hubble diagram in the form of residuals 
with respect to the specific choice of Q,m=Q-'2.1 and S1a=0.73, for 
SNe la and GRBs, together with different lines corresponding to 
different JIm, f^A pairs. It is likely that the GRB event rate fol- 
lows or increases with z faster than the global star formation rate 
(e.g., Lloyd-Ronning, Fryer & Ramirez-Ruiz 2002; Yonetoku et al. 
2003; Firmani et al. 2004), ensuring that GRBs exist up to z 
20. Even at these redshifts, GRBs are easily detectable, but the abil- 
ity of high-2; GRBs to accurately measure the Universe depends of 
course on the errors of the relevant observables. The fact that the 
bolometric fluences of GRBs with measured z do not strongly cor- 
relate with z (GGL) means that high-z bursts can be observed with 
large signal-to-noise ratios. This ensures that E^so and -Epeak can 
be measured with good accuracy also at high z's. To derive E~i, we 
need also information related to the afterglow emission. This, for 
high-z GRBs observed at a given time t after trigger, is related to 
the intrinsic emission at an earlier t/{l + z) time, when the emit- 
ted flux was stronger. It follows that the observed afterglow fluxes 
of high-z; GRBs are not much fainter than the afterglow fluxes of 
nearby GRBs observed at the same time (Lamb & Reichart 2000). 
Thus, high- and low-2: GRBs will have comparable error bars on 
the relevant quantities required to derive their luminosity distance. 
In addition, the detection of GRBs is not limited nor affected by 
dust extinction or by optical absorption by Ly-a clouds, which are 
instead important issues for high-2: SNe la. 
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